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ABSTRACT 

The X-ray transient MAXI J1836-194 is a newly-identified Galactic black hole binary 
candidate. As most X-ray transients, it was discovered at the beginning of an X-ray 
outburst. After the initial canonical X-ray hard state, the outburst evolved into a hard 
intermediate state and then went back to the hard state. The existing RATAN-600 
radio monitoring observations revealed that it was variable on a timescale of days and 
had a flat or inverted spectrum, consistent with optically thick synchrotron emission, 
possibly from a self-absorbed jet in the vicinity of the central compact object. We 
observed the transient in the hard state near the end of the X-ray outburst with the 
European VLBI Network (EVN) at 5 GHz and the Chinese VLBI Network (CVN) at 
2.3 and 8.3 GHz. The 8.3 GHz observations were carried out at a recording rate of 
2048 Mbps using the newly-developed Chinese VLBI data acquisition system (CDAS), 
twice higher than the recording rate used in the other observations. We successfully 
detected the low-declination source with a high confidence level in both observations. 
The source was unresolved (^0.5 mas), which is in agreement with an AU-scale com- 
pact jet. 

Key words: stars: individual: MAXI J1836-194 - stars: variable: other - radio 
continuum: stars - X-rays: binaries. 



1 INTRODUCTION 

MAXI J1836-1 94 is a new X-ray transient, discovered on 
2011 August 30 (|Negoro et alfeoill ). The transient was clas- 
sified as a black hole candidate, because a relativistically 
broadened iron emission lin e was revealed by Suzaku ob- 
servations (|Reis et al.l |2012| ) during its first known X-ray 
outburst. The X-ray outburst started from the canonical 
low/hard X-ray state, then it evolved to the hard interme- 
diate state and finally we nt back to the low/ha rd state. The 
outburst was classified bv lFerrigno et al.l (§012) as a "failed" 
outburst because it never entered high/soft X-ray state. 

The transient source MAXI J1836-194 was quite bright 
in radio and infrared. The Karl G. Jansky Very Large Ar- 
ray (VLA) and the RATAN-600 radio observations revealed 
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an optically-thick spectrum with a flux density of 20 - 
50 mjy between 5 and 8 GHz (|Miller- Jones et al.l l201ll : 
iTrushkin et~alll201lh . The RATAN-600 5-GHz light curve 
indicates that its emission varied frequently by a factor of 
two on the ti mescale of days. Ve ry Large Telescope (VLT) 
observations |Russell et alj 120111 ) revealed a secular mid- 
infrared brightening, reaching 57±1 mjy at 12.0 /im on 2011 
October 11, making it one of the brightest black hole candi- 
date ever detected at 10 microns. 

In view of the rare occurrence of such a bright tran- 
sient source associated with a black hole candidate, we con- 
ducted milliarcsecond-resolution VLBI (Very Long Base- 
line Interferometry) observations of the radio counterpart 
of MAXI J 1836-194 just after the mid-infrared brightening. 
In this paper, we report the first VLBI detection of its com- 
pact radio core. 
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Table 1. Summary of the VLBI observations. 



Proj. Code Date UT Range Stations Baseline Lengths Frcq. Duration Data Rate 

(yyyy-mm-dd) (hh:mm - hh:mm) (km) (GHz) (h) (Mbps) 



CHIN06A 


2011 


10- 


10 


08 


30 - 


- 13:00 


KmShUr 


1920-3249 


8.3 


4 


2048 


CHIN06B 


2011 


lD 


10 


08 


30 - 


- 13:00 


KmSh 


1920 


2.3 


4 


1024 


RSY01 


2011 


10- 


17 


Hi 


10 - 


- 18:00 


HhWbYsTrOnMcEf 


266-8042 


5.0 


2 


1024 



Table 2. The circular Gaussian model fitting results of the detected radio component in MAXI J1836-194. 



Frequency 


MJD 


TECOR 


Right 


Ascension 


Declination 


Flux Density 


Angular Size 


(GHz) 


(day) 










(mjy) 


(mas) 


8.3 


55844.450 


None 


18 h 35 


11 43? 44439 


-19° 19' 10" 4825 


6.46 ± 0.33 


^ 0.5 


8.3 


55844.450 


JPLG 


18 h 35 


11 43? 44451 


-19°19'l0"4869 


6.64 ± 0.32 


^ 0.5 


8.3 


55844.450 


CODG 


18 h 35 


11 43? 44452 


-19° 19' 10" 4868 


7.23 ± 0.30 


0.5 


5.0 


55851.712 


None 


18 h 35 


11 43? 44449 


-19°19'l0"4917 


4.09 ± 0.56 


^ 0.4 


5.0 


55851.712 


JPLG 


18 h 35 


11 43? 44451 


-19° 19' 10" 4923 


4.33 ± 0.45 


^ 0.4 


5.0 


55851.712 


CODG 


18 h 35 


11 43? 44455 


-19°19'l0"4980 


4.03 ± 0.57 


^ 0.4 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 CVN observations at 3 Gbps recording rate 

We conducted the first-epoch observations of MAXI J 1836- 
194 during a dual-frequency, wide-band test experiment 
with the Chinese VLBI Network (CVN) on 2011 Octo- 
ber 10 (project code CHIN06; Table [[]). Three anten- 
nas participated: Shanghai/ Sheshan 25m radio telescope, 
Urumqi/Nanshan 25m radio telescope, and Kunming 40m 
radio telescope (a new VLBI station in Yunnan Province, see 
lHao et al-llibioh . Each radio telescope observed only right- 
hand circular polarisation and used two back-ends: Mark4 at 
2.3 GHz and CPAS (Chine se VLBI Data Acquisition System 
System, IZhang et al]|2010l , at 8.3 GHz. Both back-ends were 
configured to their maximum recording rate: 1024 Mbps 
(16 x 16 MHz, 2-bit sampling) for Mark4, and 2048 Mbps 
(16 x 32 MHz, 2-bit sampling) for CDAS. 

We selected PMN J1832-2039 as the phase-reference 
calibrator source. Its J2000 position according to the global 
astrometric solution rfc_2010o3 is: a = 18 h 32 m ll:046479 
(la = 0.2 mas), S = -20° 39' 48'.' 20367 (la = 0.4 mas). 
It lies l.°6 away from the target, to the south-south- west. 
The duty cycle time was 3 minutes: 1 minute on the cali- 
brator and 2 minutes on the target. We observed this pair 
of sources for four hours. The remaining two hours in the 
experiment were used to observe OQ208 to verify the perfor- 
mance of the CDAS back-end. The data were correlated with 
the EVN Software Correlator at JIVE (SFXC). The corre- 
lation parameters were 1 second integration time and 256 
frequency points per subband. The Urumqi 2.3-GHz data 
were lost due to a recording failure. 

2.2 EVN observations in e-VLBI mode 

We performed the second-epoch observations with the EVN 
at 5 GHz on 2011 October 17 (EVN project code RSY01; 
see Table [1]). The participating stations were Effelsberg, 
the phased array Westerbork Synthesis Radio Telescope 
(WSRT), Onsala, Medicina, Hartebeesthoek, Yebes, and 

1 http: / / astrogeo.org/rfc 



Torun. The observing set-up was 1024 Mbps, using both 
circular polarizations. The phase-reference source remained 
PMN J1832-2039 for these observations. Furthermore, we 
also observed two nearby sources: NVSS J183542-191943 
(0.5 arcminutes away from the target, 35.1 mjy at 1.4 GHz) 
and NVSS J183522-191743 (5 arcminutes away, 28.7 mjy at 
1.4 GHz) in the hope that these could be used as secondary 
phase-reference calibrators. The duty cycle time was 5 min- 
utes: 1 minute spent on the main calibrator, 1.5 minute on 
the transient, and 1.5 minute on one of the NVSS sources 
(each alternating in successive phase-reference cycles). This 
epoch lasted for two hours. The data were streamed to JIVE 
and correlated in real-time (e-VLBI) using the EVN MklV 
Data Processor with 1 second integration time and 64 fre- 
quency points per subband/polarisation. 



2.3 Data reduction 

We calibrated the data using the NRAO Astrono mical Im- 
age P rocessing System (AIPS, version 31DEC11, iGreisenl 
2003), using the following procedure. (1) We removed the 
ionospheric delay using total electron content (TEC) mea- 
surements from GPS monitoring. (2) We calibrated the cor- 
relation amplitudes based on the measured system temper- 
atures and antenna gain curves. We used nominal system 
equivalent flux densities (SEFD) for Kunming and Shang- 
hai. (3) We corrected for the antenna parallactic angle. (4) 
We corrected the position of the Kunming telescope with 
the AIPS task CLCOR to account for improvement in its 
position estimate that became available after the correlation 
(AX = 0.157 m, AY = 0.486 m, A Z = 0.458 m). (5) We 
then fringe-fitted the PMN J1832-2039 data, and corrected 
the antenna bandpasses based on the calibr ator as well. (6) 
We i maged PMN J1832-2039 in Difmap (|Shepherd et all 
Il994f ). after all the calibration solutions were applied. (7) 
We removed the phase contribution arising from the source 
structure through a further step of global fringe fitting using 
the calibrator image in AIPS. (8) We transferred the phase 
solutions to the target sources via linear interpolation. (9) 
We self-calibrated the calibrator and applied its amplitude 
and phase solutions to the target source. Note that the am- 
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Figure 1. VLBI total intensity images of the Galactic black hole binary candidate MAXI J1836-194. Both images were produced with 
natural weighting. The synthesised beams are plotted in the bottom-left corner (left image: 3.7 X 1.3 mas at the position angle 21.8 
degree, right image: 6.1 X 4.6 mas at the position angle —12.9 degree). The contours start from 3<r off-source noise level (0.54 mjy beam - 1 
at 8.3 GHz, 0.24 mjy beam -1 at 5 GHz) and increase by a factor of 2. 



plitude self-calibration was not possible in the 8.3 GHz ex- 
periment because of the small number (< 4) of the stations. 
We then averaged the data in frequency over all subbands. 
(10) We identified the transient source in the dirty map and 
measured its position by fitting a point source model to the 
visibility in Difmap. (11) After the model position was fixed, 
we self-calibrated the data using a 2-minute solution inter- 
val to improve the noise level and look for possible other 
components. We kept Hartebeesthoek (long baselines) and 
Onsala (less sensitive) phases fixed in this process. 

We also reduced the 5-GHz WSRT synthesis array data 
in AIPS, using 3C 286 to calibrate the flux densities. We 
self-calibrated PMN J1832-2039 (an ideal point source for 
the WSRT) in both phase and amplitude, and applied the 
solutions to the target. Finally, we averaged the data in each 
subband and imaged the source in Difmap. 



3 VLBI IMAGING RESULTS 

The VLBI imaging results are displayed in Figure [T] 
MAXI J1836-194 is unresolved at both frequencies; the visi- 
bility data can be fitted with a point source model very well. 
The model-fitting results are listed in Table [2] The phase- 
reference calibrator image is shown in Figure [2] The cali- 
brator source had a bright, one-sided jet. The eastern com- 
ponent has the most compact structure with similar peak 
brightness (~0.38 Jybeam -1 ) at the two frequencies. Thus, 
we assume this to be the radio core. We measured the source 
position with respect to the compact core of the calibrator 
at each frequency. Note that this new position of the target 
sourc e is different from the one we reported earlier l|Xu et al.l 
l201ll ): in that telegram, the reference origin was the central 
bright jet component rather than the radio core. In order to 
estimate the upper limit of the angular size of MAXI J1836- 
194, we also fitted the visibility data with a circular Gaussian 
model. The last column in Table [2] presents the estimates of 
the angular size. Both observations give a consistent upper 
limit of 0.5 mas. 

We also searched for extended emission in the VLBI im- 
age by increasing the weights of short-baseline data. There 



was no extended emission found. Based on the WSRT obser- 
vations, the total flux density of the transient at 5 GHz was 
5.7 ± 0.2 mjy. Therefore we recovered 95% of the total flux 
density in the simultaneous VLBI observations. If there is 
extended emission, its contribution to the total flux density 
is quite low, ^ 5%. 

There was only one baseline (Kunming - Shanghai) at 
2.3 GHz. The baseline gave a total flux density 5±1 mjy. Be- 
cause phase self-calibration was not possible, there may have 
been a significant loss of phase coherence at this frequency 
(note the very low declination). Assuming ~30% loss, similar 
to what we found after applying the phase self-calibration 
at 8.3 GHz, the source would have a flux density close to 
that at 8.3 GHz. Thus, the source did not have a spectral 
index significantly different from on the date of these dual- 
frequency observations. 

The two nearby NVSS sources that we included in 
the 5 GHz observations as candidate secondary calibra- 
tors within 5' of MAXIJ1836-194 were not detected at the 
threshold of 5a ~ 1.3 mjy. Therefore, they could not be 
used to estimate spatial gradients of propagation effects to 
the phase referencing in the vicinity of the target. 



4 DISCUSSION 

4.1 MAXI J1836-194 position stability 

After careful removal of source-structure effects in the phase- 
referencing (c.f. Sect. 12.31 steps 6-8), there is still a notable 
position difference (Aa=0-1.4 mas and A5=5.4-11.2 mas) 
between 5 and 8.3 GHz. There may well be a frequency- 
dependent core shift in the reference source J1832-2039 that 
could not calibrate absolutely with available observations, 
but this would lie along the E-W jet axis (see Fig. [2| and 
generally be on the order of 0.1 mas between 5-8.3 GHz 
jKovalev et ai1 l2008). Such a core shift would not explain 
the apparent position difference which is predominantly in 
declination. Other critical characteristics of these observa- 
tions include the short duration and the very southerly tar- 
get (and an even more southerly reference source), leading to 
persistently low elevations from European VLBI antennas. 
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At these low elevations, even the l.°6 separation between the 
target and reference sources is larger than ideal. Thus the 
first suspects would be imperfectly calibrated propagation 
effects in the phase-referencing. 

These observations from October 20fl took place at a 
time of relatively elevated solar activity coming out of the 
prolonged quiet period, which could lead to a more active 
ionosphere. Quantitatively, the daily solar 10.7cm flux den- 
sity (e.g., NOAA National Geophysical Data Centre), aver- 
aged over a 27-day solar rotation period was about 140 [10 4 
Jy], a value not seen since late 2003. The EVN MklV Data 
Processor applies no a priori ionospheric delays; rather it 
is standard to apply IONEX maps f r om the International 
GNSS Service (IGS) (|Dow et alj|2009l : I Walker fc Chatteried 
ll999T l in AIPS data reduction, as described in Section 12.31 
There are a few IGS analysis centres who process obser- 
vations of GPS satellites from 100-150 ground stations to 
produce these IONEX maps of vertical column density of 
electrons, gridded spatially at (A$ = 2?5, AA = 5?0) and 
temporally at 2-hour intervals. The AIPS task TECOR com- 
putes phase corrections for the lines of sight from the VLBI 
antennas to the sources throughout the observations by in- 
terpolating within these IONEX maps. Table[2]shows the ef- 
fect of treating the ionospheric delays in three ways: making 
no correction ("None"), and then using the IONEX maps 
from the analysis centres at JPL ("JPLG") and the Uni- 
versity of Bern ("COGD"). Using IONEX maps from these 
different Analysis Centres leads to changes of up to 6 mas 
in the declination of the target, which is on the order of the 
observed position difference. 

The correlator model computed dry and wet tropo- 
spheric d elays from the Saa stamoinen zenith-delay formula- 
tion (e.g.. lDavis et al.lll985l ). parameterized by surface pres- 
sure, temperature, an d relative h umidity, and the Neill map- 
ping function (NMF. [Neilllll996h to convert from zenith to 
slant delays along the lines of sight. At low elevations, mis- 
modelling of the zenith delays can become amplified (sec z 
for the RSY01 observati on average 2-5 fo r the v arious Eu- 
ropean VLBI antennas b iBrunthaler et alJ l|2005h developed 
a method to determine improved zenith-delay estimates via 
scheduling geodetic- VLBI-like blocks into the observations, 
but the short duration of the present experiments precluded 
use of this tactic. To characterise the range of position shifts 
that may be induced by mismodelled tropospheric delays, 
we computed zenith delays in two ways: using the surface- 
weather values for antennas that had such informatio n in 
their log files and using Calc-10 (|Gordon et alJ I2006T I de- 
faults. We applied the NMF to both, and computed differ- 
ential tropospheric delays. The median of these was up to 
~50 ps, corresponding to a phase offset ~100° at 5 GHz. 
Through adjusting the position of the target and running 
phase self-calibration, we found that similar phase change 
results from a position shift of ~5 mas. 

Thus the variation in the position of the target esti- 
mated from the CHIN06 and RSY01 observations is con- 
sistent with effects of ionospheric/tropospheric propaga- 
tion in the phase referencing. Similar effects were also no- 
ticed in other low-declin ation sources, such as H 1743-322 
l|Miller- Jones et al J 12012T b 

Nonetheless, we can also briefly consider physical causes 
of a position shift associated with the target itself over the 
week between the two observations. The transient produced 



PMN J1832-2039 EVN at 5 GHz 
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f " % f g i Q Total Flux Density: 0.99+/-0.1 Jy 
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Figure 2. VLBI total intensity images of the phase-referencing 
source PMN J1832-2039. Both images were produced with natu- 
ral weighting. The synthesised beams are plotted in the bottom- 
left corner. The contours start from 3<r off-source noise level (0.5 
mjy beam - 1 at 5 GHz, 1.0 mjy beam - 1 at 8.3 GHz) and increase 
by a factor of 2. 

a "failed" outburst (|Ferrigno et al.ll2012T ). never transiting 
to the high/soft X-ray state during its evolution (low/hard, 
hard intermediate, and back to low/hard X-ray states). Both 
of our observations were carried in the canon ical low/hard 
state , where major ejection eve nts are unlikely l|Corbel et al.l 
|2004 iFender et al.ll2004l2009l) . Moreover, our target source 
was quite compact at both epochs with brightness tem- 
peratures ~10 8 K, showing no sign of extension between 
the two epochs. That the observed spectrum between 2.3 
and 8.3 GHz is consistent with being optically thick indi- 
cates that we detected a compact, synchrotron self-absorbed 
component of the system — its radio core. This is in agree- 
ment with the fast va riability observed by the RATAN-600 
l|Trushkin et alJ 1201 ll ) during the hard state at an earlier 
epoch. Thus the position shift is unlikely to be related to 
moving ejecta. 

The position measured at 8.3 GHz is likely more reli- 
able on several grounds: the beam pattern of the 8.3 GHz 
data was clearly seen in the dirty map, the underlying CVN 
observations had twice the on-source time and elevations 
that were generally 10° higher, and the higher frequency 
affords 40% less ionospheric-phase contribution for compa- 
rable TEC levels. 

4.2 VLBI detections as the compact-jet region 

The arguments presented above convince us that the point- 
like detections in Figure [T] correspond to the same object, 
i.e., the radio core of MAXI J1836-194. The existence of a 
synchrotron jet in the vicinity of the stellar-mass compact 
object in MAXI J 1836-194 is furt her supported by the opti- 
cally thick spectrum in the radio (jTrushkin et al.|[201ll ) , and 
the optically thin p ower-law spectrum in the mid-infrared 
l|Russell et al.l 1201 lT ) . Similar compact jets have been di- 
rectly resolved in t he hard state in blac k hole X-ray binaries, 
e.g., Cygnus X-l ^Stirling et all 1200 ll ) and GRS 1915+105 
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(|Dhawan et all l2000h . In the present case of MAXI J1836- 
194, the compact jet remains unresolved. We can constrain 
its angular size to be ^0.5 mas, which corresponds to an AU- 
scale linear size for a Galactic source. Because MAXI J1836- 
194 lies outside of the Galactic plane (b = —5? 4), this con- 
straint should not be affected by scatter-broadening arising 
from the Galactic medium along the light of sight, which 
could potentially smear out the intrinsic source structure at 
low frequencies. We compare the apparent size constraint 
with the known, resolved compact jets below. 

Within the fr amew ork of the jet model of 
iFalcke fe Biermannl (|l999T ) with minimal number of 
free parameters, the characteristic size #j e t of the compact 
jet significantly depends on the mechanical power of the 
jet: djet oc Qj' c f xi , where Xi ls a parameter depending on 
the inclination angle of the jet, with a value generally 
around 1, and the jet power Qj e t is on the order of the 
X- ray luminosity of the coupl ed disk-jet system considered 
by IFalcke fc Biermannl (| 19991 ). MAXI J1836-194 had an 
X-ray luminosity ~ 10 36 ergs -1 between .3 - 10 keV 
(corrected for absorption, iKennea et al.l 1201 lh assuming a 
distance of 8.5 kpc, three orders of magnitude lower than 
that (~ 10 39 ergs -1 ) of GRS 1915+105. Assuming that 
the parameters not determined from these observations 
(distance, mass and inclination angl e) are the same as 
GRS 1915+105 (|Dhawan et all l2000h . MAXI J1836-194 
would have a characteristic size at least an order of mag- 
nitude smaller than that of GRS 1915+105 (~0.7 mas at 
8.3 GHz). This result is consistent with our upper limit 
of the source size. We therefore conclude that we have 
detected the unresolved, synchrotron self-absorbed jet base 
in MAXI J1836-194, in accordance with the expectations 
for its hard/low X-ray spectral state. 



5 CONCLUSIONS 

We have presented VLBI observations of the recently de- 
tected Galactic variable X-ray source MAXI 1836-194 with 
the EVN at 5 GHz and the CVN at 2.3/8.4 in October 2011. 
Despite the unfavourable celestial position for the telescopes 
located in the Northern hemisphere, this low-declination 
Southern object has been firmly detected at all three fre- 
quencies as an unresolved radio source. The VLBI data allow 
us to put a firm upper limit on the angular size of the source 
as 0.5 mas, which corresponds to an AU-scale linear size for 
a Galactic object and is as expected within the framework 
of compact synchrotron jet models. 
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